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Abstract 

We study the Lie point symmetries of Einstein's equations for the 
Friedmann-Roberstson-Walker Cosmology. They form either a two - di- 
mensional or a three - dimensional solvable group depending on the form 
of the self interacting potential. Using the invariants of the group we re- 
duce the second order system of differential equations into a first order 
system. Writing the action in terms of the proper time we study the point 
symmetries and the variational symmetries of the resulting equations. 
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1 Introduction 

A symmetry group of a system of differential equations is a Lie group acting on 
the space of independent and dependent variables in such a way, that solutions 
are mapped into other solutions. Knowing the symmetry group allows one to 
determine some special types of solutions that are invariant under a subgroup 
of the full symmetry group, and in some cases one can solve the equations 
completely. The symmetry approach to solving differential equations can be 
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found, for example, in the books of Olver [10], Bluman and Cole [1], Bluman 
and Kumei [2], Fushchich and Nikitin [5] and Ovsiannikov [11]. 

Einstein's General Theory of Relativity is based in the most fundamental 
way on the concepts of symmetry and general covariance. Lie's Theory on the 
other hand is the most systematic mathematical way to study symmetries of dif- 
ferential equations. So it is obviously interesting to apply the Lie methods to a 
theory fundamentally related with the idea of symmetry and general covariance. 
In this paper we study the Lie point symmetries of the notorious Einstein's equa- 
tions for the Friedmann-Roberstson- Walker Cosmological model with a matter 
field of self interacting potential V((f>). For completeness wc include here their 
derivation from the Einstcin-Hilbert action. For the above Cosmological model 
it is given by 

1 f , N T f a da\ 2 , 2 { a 2 d<p\ 2 „ 4 i 

S = 77 dt— - — -r +ka 2 +[—-i-) - 2a*V 1 

2 J a L \N dt) \N dt) J K ' 

where N is the lapse function and a is the scale factor or the size of the Universe, 
and (j) is the scalar field of self interacting potential V{4>). If we vary the action 
S with respect to N, <j> and a we will get the following equations: 

a 2 = 2a 2 V + a 2 ft - k (2) 

= + 3^+^ (3) 
a dep 

Q = -2ad-a 2 -k-Sa 2 ft + 6a 2 V (4) 

where the dot above means ^ . Since we have imposed the Hamiltonian con- 
straint we can set N = 1. Using Equations (2) and (3), Equation (4) becomes: 

a = 2aV - 2a(j) 2 (5) 

which is usually called the second Einstein equation. Here k — +1 for the case 
of a closed universe for which the space part is a three sphere, k = for the case 
of a marginally open universe with flat space part and finally k = — 1 represents 
an open universe whose space part is a three dimensional hyperboloid embedded 
in Minkowski space. 

The above equations are not only important for classical cosmology but they 
are also important for the semiclassical approximation to the Wheeler-DeWitt 
equation [4], i.e., to the semiclassical quantum cosmology. The Wheeler-DeWitt 
equation for the Friedmann-Roberstson- Walker Cosmology takes the form: 

(-iv 2 + [/)* = (6) 

where the Laplacian operator is calculated in the minisuperspace metric 

ds 2 = -ado 2 + a z d(j) 2 (7) 
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and the potential is given by 



U = - a -+a*V (8) 

If we write the wavefunction of the universe in the form 

* ~ Ce iS (9) 

where C is slowly varying with respect to 5*, then in the W.K.B. approximation 
we choose the action to satisfy the Hamilton- Jacobi equation 

I(VS) 2 -^+a'V = (10) 

which is consistent with the fact that '5 satisfies the Wheeler-DeWitt equation. 
The action S defines then a vector field 

| = VS.V (11) 

It can be easily proved that its trajectories satisfy the above classical equations 
(2), (3) and (5). See for example [12]. In a forthcoming paper we study the 
symmetries of the Wheeler-DeWitt and of the Hamilton- Jacobi equations and 
we compare them with the results of this work. Finally we note that the quantity 
a 2 — a 2 <j) 2 — 2a 2 V(4>) is a constant of motion for the second order equations (3) 
and (5), which we study next. 



2 Study of Lie point symmetries 

We study here the Lie point symmetries of the equations 

a = 2aV{(j)) - 2a0 2 (12) 

* = d3) 
a dip 

where V((f>) is an arbitrary function of <p. Here we give some details of the 
calculations, because later we apply the same method to equations (3) and (4) 
without showing any calculations. If the vector field G is given by 

d d d 

then its second prolongation or extension is given by 

d , d r d . ., d , i . • . d 



pr^'G =t— + A— + + (A - ra — + ($ - r^-y 
ot da dip da dip 

+ (A - fa - 2fa)— + ($ - f (p - 2i'ij)) — 
da dij) 



(15) 
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In the above equations t, A and <I> are functions of t, a and <fi. The necessary 

and sufficient condition for G to be an infinitesimal generator of the symmetry 
group of the above equations (12) and (13) are: 

pr {2) G(Eql2) = (16) 

pr (2) G(Eql3) = (17) 

or equivalently 

A-fa-2fa = 2AV + 2A$^ - 2A/> 2 - 4a0($ - f </>) (18) 



(19) 



(20) 



$ - f4> - 2f'<j> = -3(i -fa)-- 3(4 - f(£)- 

Expanding Equation (18) we get 

A tt + 2A ta a + 2A t4> <f> + 2A a4> aj> + A aa a 2 + A^ 2 + 2aV{(p)A a 

- 2a4> 2 A a - 3—^ - - AaV{4>)T t - AaV{4>)T a a 

- AaV{4>)T 4> 4> + Aa4> 2 T t + Aa4> 2 T a a + 4a<fi 3 T$ - T tt a - 2r ta d 2 

- 2r t ^(j)a - 2T a<j> a 2 4) - r aa d 3 - T^cj) 2 a - T a a2aV((j)) + T a a2acj) 2 

+ 3r ^ + r^%- = 2AV(<j>) + 2a$^ _ 2 A^ 
a d<p d<p 

- 4a^$ t - Aacf)<& a a - 4a^ 2 $0 + 4a0 2 r t + 4ar a (j) 2 a + iar^ 3 

Here we follow the same method as in [9] . On equating the coefficients of d we 
get r aa = 0, thus 

T(t,a,(f>) =c l (t,4>)+c 2 (t,4>)a (21) 
From the terms of a 2 <f) we obtain 

- 2r a(jl + 3^ = (22) 

which using Equation (21) and after some algebra gives 

t = ci(t) +c 2 (t)a (23) 

The terms of atp 2 give 2ar a — t^, thus from the previous equation we infer 
that r a = 0, so t depends only on t. Checking the coefficients of d 2 we get 

A = c 3 (t,ct)) + c 4 {t,(j))a (24) 
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From the terms of <j> 2 and using Equation (24) we obtain: 

w +2C3+ w a+4a ** =0 (25) 

Similarly from the coefficients of a<j> and using again Equation (24), after some 
algebra we end up with the following relation: 

On differentiating Equation (25) with respect to a we get 

d 2 CA 

a ~d^ + 4a * + 4a2 * a< ^ = ( 2? ) 
Similarly from equation (26) by differentiation with respect to <f> we have: 

Equations (27) and (28) give: 

2a w +3 w a+Aa ^=° (29) 

which combined again with Equation (25) implies 

<9 2 c 4 „<9 2 c 3 „ , . 

a w = ~ 2 w +2cs m 

This last equation gives the functions c%(t, cf>) and Ci(t, 4>) in the form: 

C3(t,0) = h(t)e + fc(t)e^ (31) 

c 4 (M)=/(i) + <K^ (32) 

where at the moment ft, fc, / and 5 are arbitrary functions of time. Substitut- 
ing (31) and (32) back in Equation (25) we get for the function $>(t,a,4>) the 
following result 

* = -f«' + f«- + *,<) (33) 
Finally Equation (26) combined with (31), (32) and (33) gives: 

p(a,t) = ^g{t)lna + iJ, (34) 

Here we put together what we have found up to now and we get the expressions 
below for the functions A(t, a, <f)), $(t, a, </>) and r(i, a, <f>): 

A = he 4 ' + he- 4, + [/(*) + g<j>]a (35) 
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9 = ~he* + ^-ke-* + \g]n.a + fi (36) 
4a 4a 4 

t = r(f) (37) 

Next from the coefficients of d, <ft and after some algebra we end up with the 
relations: 

g=h=k= ft (38) 

and 

2/ = f (39) 
The terms that do not involve derivatives of a or (f> give 

A a + 2aV((f>)A a - - AaV{4>)f = 2AV(<P) + 2a$^ (40) 

Substituting (35), (36) and (37) into the last relation and looking at the terms 
that do not depend on a after some algebra we get h = k = 0. Similarly, the 
term proportional to a In a gives 5 = and from the terms linear in a we end 
up with the relation 

dV 

f - AV(4>)t = 2—ii (41) 

From the last relation we infer that, if ^ = 0, or equivalently if V((j>) is pro- 
portional to e -2 * 1 we must have f = ji and / = 0, from where we get r = fj,t + c 
and /(f) = cif + c 2 . Then using Equations (38) and (39) we get c\ — 0. Putting 
all these together the above equations (35) through (37) become: 

A = aa (42) 

$ = fi (43) 

r = /if + c 2 (44) 

when V{4>) = ce~ 2< ^ . On the other hand if V{4>) ^ ce~ 2 * Equation (23) gives 
/ = f = = and using equations (38) and (39) we get / = 0. So we have 

A = cia (45) 

$ = (46) 

r = c 2 (47) 

It is straight forward to prove that the vector fields determined by the Equa- 
tions (42)-(44) and (45)-(47) both satisfy the condition just mentioned, i.e., they 
arc symmetries of the Klein-Gordon equation. So we have the following theorem. 

Theorem 1. The Lie point symmetries of the ordinary differential Equations 
(12) and (13) with V(<p) = ce -2 ^ is a three dimensional Lie group with generator 
G = (M* + c 2 )i+cia^+ M ^ 

On the other hand if V(<f>) arbitrary and ^ ce~ 2 ^ we get a two dimensional 
Lie group generated by 

c 2i + c i a lr 

which is a subgroup of the previous one. 
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It is interesting to note that in the case where the self interacting poten- 
tial V(4>) is exponentially decaying the Hartle-Hawking boundary conditions 
in quantum cosmology take an especially simple form [7]. We do not know 
if there is any connection between this result and the above theorem , where 
again the exponentially decaying potential plays a special role. The quantity 
E = a 2 — 2a 2 V((f>) — a 2 <j) 2 , which is — k from the Einstein's constraint equation, 
is preserved by the solutions of the above Equations (12) and (13). We can 
easily find how the above vector fields act on E: 

pr^G(E) = 2(c 1 -iJ l )E (48) 

This shows that in the case of a flat universe G is a symmetry group for all 
Einstci's Equations (2), (3) and (5). The generators of the symmetry group G 
are given by: 

Z = al (51) 

The vector field X corresponds to c\ = c 2 = 0, /i ^ 0, the vector field Y 
corresponds to c\ = /i = 0, c 2 7^ and Z corresponds to c-i = /j, = 0, c\ ^ 0. 
Their multiplication table is 

[X,Y] = -Y (52) 

[X,Z]=0 (53) 

[Y, Z] = (54) 

So here we have a solvable but not nilpotent Lie group. For closed universes 
and for open hyperbolic only, Y is a symmetry for all Einstein's equations. In 
the special case of V((f>) = ce~ 2< ^ the subgroup generated by Y and W = X + Y 
leaves invariant all Einstein's equations for all the cases where k = — 1, 0, +1. 

In a recent study of the Maxwell-Bloch system [3] , we found again only solv- 
able Lie groups of symmetry. We do not know if there is any general condition 
which both systems satisfy and which forces them to accept only solvable groups 
of symmetry. This has to be investigated. 



3 Reduction 

It is well known that if a one parameter Lie group of transformations is ad- 
mitted by an ordinary differential equation then its order can be reduced by 
one. For first order ordinary differential equations this corresponds to a re- 
duction to quadrature. This reduction of order can always be accomplished by 
using canonical coordinates accosiated with the group. For higher order ordi- 
nary differential equations the reduction in order can be accomplished by using 
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differential invariants. Here we will reduce the second order system consisting of 
the Equations (12) and (13) to a first order system using the change of variables: 



<j> = x (55) 

4> = V (56) 

= w (57) 

which are invariants of the symmetry group generated by Y and Z. It is trivial 
to show that under the above transformation the second order system (12)-(13) 
reduces to the following first order system of ordinary differential equations 

d y o dV (zo\ 

y Tx = ~ 3wy ~ to (58) 

y ^+ w 2 =2V(x)-2y 2 (59) 
ax 

and two quadratures. If we find y{x) and w(x) then we can find <f>(t) and a{t) 
by quadratures: 

/' ,/, , 

t (60) 



/ 



y{4>) 

Inverting the relation above we can get as a function of t. Then we can find 
a as a function of t using the quadrature 

a = cxp J w((f>(t))dt (61) 

Unfortunately the reduced system consisting of the Equations (59) and (60) is 
not autonomous even though the original one is. Using x, y and w the quantity 
E can be written in the form 



E = a 



w 2 - 2V(x) - y 2 (62) 



and it is trivial to prove that the total derivative of E with respect to x is zero. 
So E is still conserved even though is not expressed solely in terms of the new 
variables x, y and w. Both Equations (12) and (13) are locally solvable and of 
maximal rank since their rank is the rank of the matrix 

/ -2V + 24> 2 -2a^ 4a4> 1 \ 
I "3# 3* l) 

As a result, conditions (16) and (17) are necessary and sufficient for G to be a 
Lie point symmetry. 
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4 Variational symmetries 



Even though Equations (2), (3) and (5) can be reduced from a variational prin- 
ciple, as we have seen in the introduction, the study of variational symmetries is 
problematic and the application of Noether's theory to this specific variational 
process, to our knowledge, is an open problem. 

We will look here more closely at the peculiarities of our variational prob- 
lem. As we have seen in the introduction, the original variational equations are 
Equations (2), (3) and (4). We have replaced Equation (4) by Equation (5) 
since combining algebraically Equations (2) and (4) we get (5). Thus we can 
use Equations (2), (3) and (5), which are equivalent with the original equations. 
These are three equations for the three uknown functions N(i), a(t) and <j>(t). 
But actually they are not independent, because using (2) and (3) we can prove 
(5); using (2) and (5) we can prove (3) and even using (3) and (5) we can prove 
(2) (but we can not infer the value of k). On the other hand if in the above 
equations we rescale the time, by using as our new time the integral of the lapse 
function, i.e., the proper time of general relativity, we will have three (depen- 
dent) equations for only two unknown functions, a(t) and <t>{t), where t now 
stands for the proper time. The resulting equations are not in any obvious way 
the Euler-Lagrange equations of a variational problem and from here it springs 
the difficulty of applying Noether's theory. It looks that N is a quantity we 
need to formulate the problem but eventually disappears from the dynamics, it 
is like an ignorable quantity. This is related with the problem of distinguishing 
the gauge and the dynamics in general relativity a main source of difficulties for 
both classical and quantum general relativity. See for example the articles of 
Teitelboim [13], Kushar and [8], Hanson, Rcgge and Teitelboim [6]. 

Here we will use from the beginning the proper time in the Einstein-Hilbert 
action and we will study the resulting variational problem. The action takes 
the form 



S=l I dt 



^) 2 + ka + a 3(^)-- 2a z V 
dt J V dt 



(63) 



where t from now on will represent the proper time. Varying with respect to 
<p and a we get the Equations (3) and (4) of the introduction. We do not 
have here the Hamiltonian constrained equation and as a result we do not have 
also Equation (5). In short the resulting equations are not equivalent with the 
system we had before. Following the same method as in Section 2 we get the 
following result: 

Theorem 2. The Lie point symmetries of the Equations (3) and (4) with k = 
is a two dimensional solvable Lie group with generators G = C2 + c \ a ~§^ For 
k ^ the only symmetry is time translation. The above hold independently of 
the form of the potential V(<p). 

Equations (3) and (4) are the Euler-Lagrange equations for the Lagrangian 
L = — iad 2 + ifca + ia 3 </> 2 — a 3 V We can easily prove that the time translation 

is a variational symmetry, since it satisfies the condition pr^ 1S) X{L) + LSI ■ £ = 
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which is necessary and sufficient condition for a vector field X to be a variational 

symmetry. On the other hand scaling in a is not a variational symmetry. From 

Noether's theory we can find the conserved quantity associated with time trans- 

_ v i 
lation: If X = £ A* J,+ £ 1> a (t 

' u ) du a ^ ^ variational symmetry, where t* are 

i— 1 a— 1 

the independent and u a are the dependent variables. Then the characteristics 



of X are given by Q a = ij: a 



£ A 1 and they generate a conservation law in 



characteristic form by the equation V • P = £ Q a ■ E a (L) For the variational 

a=l 

symmetry Jj we can easily find its characteristics Q a — —a and — —<j>. The 
conservation law takes the form ^ = Q a ■ E a (L) + • E^{L). where L is the 
above Lagrangian and E a , E^ are the Euler-Lagrange operators corresponding 
to a and <f>. Substituting from Equation (63) we get: 



Qa ■ E a (L) + ■ E^(L) 



2~dl 



a(a 2 + ka- 2a i V 



,3 12 



(64) 



So P = a ■ E is a conserved quantity for Equations (3) and (4), where E is 
the conserved quantity we have for Equations (12) and (13). We could find the 
above conserved quantity P by considering t as a function of cither a or <j). If a 
is our independent variable the action takes the form 



dat a 



a 

*2 



ka 



2a A V 



(65) 



the Lagrangian does not depend on time, as we would expect, and the conjugate 



momentum is conserved. It is trivial to check that this conserved quantity 
agrees with P. Actually this is the Hamiltonian related to the Lagrangian above, 
something we should expect since it is related to the time translation symmetry. 

For k = we get the same symmetry group as in sections two and three, so 
we can use as new variables the invariants x, y and z of section three to reduce 
the order of the Equations (3) and (4). Indeed, in terms of x, y and z they take 
the form 

dy dV 
y— — h 3wx H — — 



dx 



dx 







2y 



dw 
dx 



3w 2 + 3j/ 2 - 6V(x) = 



5 Variational symmetries for the full system 

In the original action the Lagrangian is written in terms of a, <j), N and t, which 
is the affine time, in the form 

L = -\^ 2 + \ kNa + Na " v{<t>) (66) 
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The equations of motion for a, <f> and N take now the form 

a 2 + N 2 k - a 2 ^ 2 - 2N 2 a 2 V = (67) 

aN4> + 3Na<j> - aNj) + N 3 a^- = (68) 

d<p 

and 

Nh + 2Na(j) 2 -Na- 2N 3 aV = 0) (69) 

We can easily verify that the vectors X, Y, Z of section three are still Lie point 
symmetries for Equations (66), (67) and (68), i.e., they satisfy the following 
conditions 

pr ( 2)X(Eql) = 
for k = and V{(f>) — cexp(— 2(f>) and 

pr^X(Eq2) = pr^X(EqS) = 
for arbitrary V((f>) and k. Also, 

pr {2) Y(Eql,Eq2,Eq3) = 
for any V(<f>) and k. And finally, 

pr^Z(Eql) = 

for any V(<f>) but k = and 

pr {2) Z{Eq2 1 Eqi) = 
for any V(</>) and k, where now 

pr^X = td t +8^- ada - (pdj, - Nd^ - 2ad a - 2^ - 2Nd& 

pr^Y = d t 

pr {2) Z = ad a + ad a + add 

Here we can also check which of the above vector fields are variational symme- 
tries. It turns out that X and Z are not variational symmetries, but Y is a 
variational symmetry satisfying the condition 

pr (1) Y(L) +LV = 

The characteristics of Y are Q a — —a, — —<p, Qn — —N and applying again 
Noether's theorem we get the corresponding conservation law in the form 

^ = QaEa(L) + Q+E+{L) + Q N E N (L) 

which gives 

^ = ± \ (a 2 + N 2 k - 2N 2 a 2 V - a 2 ^ 2 ] 
dt dt 12N\ Y Jl 

This agrees with the conserved quantity P of section four if we set N = 1 . Of 

course K does not give a new conservation law since using Equation (66) we get 

K = 0. 
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